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Abstract. The BATSE mission aboard CGRO monitored the whole sky in the 20 keV - 1 MeV energy band 
continuously from April 1991 until June 2000. Although BATSE had very poor intrinsic angular resolution, the 
data can be used to survey the entire soft gamma-ray sky with < 1° angular resolution and ~mCrab sensitivity by 
using the Earth occultation method. This method determines flux by measuring the step in the count rate profile 
in each BATSE detector as a source rises above or sets below the Earth's limb. A maximum likelihood imaging 
technique can then be used to build up all-sky maps from the images of the Earth's limb produced by occulting 
sources. However, since the Earth seen from BATSE has a radius of ~70°, the limb images that intersect at the 
positions of bright point sources have a significant effect over the area of the all-sky map. A method for performing 
image cleaning on likelihood data has also been developed and is used to effectively remove artefacts from the 
all-sky maps. This paper describes the 'LIMBO' imaging technique and presents preliminary all-sky maps of 25 - 
160 keV emission, the first to be made since the HEA01-A4 mission of 1978-79. 
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1. Introduction 

Conventional imaging techniques, of focusing light with 
^ reflecting or refracting surfaces, are not possible in the 
^> , hard X-ray or gamma-ray regions of the electromagnetic 
^ spectrum due to the highly penetrating nature of these 
5^ photons (however, see lvon Ballmoos et al.l l)200lh for a re- 
port on the first use of a Bragg diffraction gamma-ray 
lens in astrophysics). For photons above ~ 10 keV, where 
grazing incidence techniques become highly inefficient, it 
is necessary to make images by transforming the spatial 
information of the object onto the detector plane with 
some sort of modulation or coding technique that can be 
used to form a reconstructed image. 

The ability to form an image has many advantages over 
simple measurements of source fluxes made with point- 
ings of basic coUimated instruments. The ability to locate 
a source with an image allows comparison with measure- 
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ments made at different times, with different instruments 
at different wavelengths. Crowded fields can be resolved, 
minimising the possibilty of source confusion and the na- 
ture of extended or diffuse emission, such as that which 
may be seen from Supernova Remnants (SNR) or the 
Galactic Plane (CP), can be investigated. Imaging detec- 
tors can also have a large gain in sensitivity over simpler 
instruments by reducing the amount of time needed to 
make accurate observations; It is possible to make an es- 
timation of the background in the region of a source at 
the same time as the observation, negating the need to 
perform control observations away from the source. In an 
instrument with a large Field of View (FOV) it is possible 
to observe many sources in a single observation, which is 
often referred to as the multiplex advantage. 

The natural progression of the last point above is to 
image the entire celestial sphere and perform an all-sky 
survey. Knowing the locations and properties of a large 
number of sources is important in all energy bands of as- 
tronomy if only for the purposes of design, planning and 
proposal justification of other instruments and observa- 
tions. Even a survey with crude angular or spectral resolu- 
tion can be highly important; as long as the sources in the 
survey catalogue have well determined uncertainties then 
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Fig. 1. The Compton Gamma- Ray Observatory and one of the 8 BATSE detector modules (drawings courtesy of the 
National Space Science and Technology Centre, Huntsville, AL, USA). 



other observations by higher resolution instruments are 
possible. Contemporaneous measurements of a large num- 
ber of sources are important in c onducting mor phological 
studies of famihes of objects fe.g. lWestmore et al.l (|2000j)). 
This is possible for an instrument with a large FOV, which 
is able to conduct a survey of such sources free from the 
biases that may be introduced by conducting many obser- 
vations on different smaller FOV instruments at different 
times. Similarly, by observing a large part of the sky, it is 
possible to maximise the probabilty of detecting unknown 
or transient objects, such as Gamma-Ray Bursts (GRBs) 
or X-ray novae, in a fixed time. The nature of hard X-ray 
and gamma-ray imaging systems means that they often 
have much larger FOVs than is possible at optical wave- 
lengths. In the extreme case it is possible to observe the 
entire celestial sphere at once. 

In high energy astronomy all-sky im ages have been 
provi ded between 0.1-2 keV by R OSAT llSnowden et"all 
1997tl and > 1 MeV by COMPTEL llOberlack et al.lll99(it 
Schonfelder et al.l200(ll) . but there is sparse information in 
the intervening band, provided almost exclusively by the 

IIEA01-A4 experiment of 1978-79 (Levine et al 1984) . 

This experiment consited of two Nal/Csl phoswich detec- 
tors arranged so that they observed simultaneously in op- 
posite directions and spun about an axis perpendicular to 
the optic axis and in the plane of the ecliptic. This meant 
that it took 6 months for the whole sky to be observed 
and so most areas of the sky were only sampled two or 
three times during the 16 month mission. The IIEA01-A4 
survey catalogue contains 72 sources in four broad energy 
channels ranging from 13 - 180 keV, at a flux sensitivity of 
14 mCrab. The A4 experiment suffered from some prob- 
lems with changing efficiencies between the two detectors 
and the sensitivity in the high energy band was limited; 
only 14 of the sources were detected at energies greater 
than 80 keV. Although the IIEA01-A4 was an important 
mission the usefulness of the survey is limited by the low 



spectral range and resolution and the fact that many tran- 
sient sources were missed because of the slow sampling of 
the whole sky. 

The aim of this work is to show how the wealth of the 
BATSE data set, which comprises 9 uninterrupted years 
of all-sky monitoring, can be exploited to achieve the max- 
imum scientific gain. 

2. The Burst and Transient Source Experiment 

The Burst and Transient Source Experiment (BATSE) 
monitored the whole sky continuously from April 1991 - 
June 2000 as an experiment aboard the Compton Gam ma- 
Ray Observatory (CGRO) (,Fishman et alJ Il989a|) . It 
was designed primar ily to detect and locate GRBs 
(Fishman et~aI] ll989bV BATSE consisted of eight uncol- 
limated 2025 cm^ Nal(Tl) Large Area Detectors (LADs), 
which had a combined field of view of An steradians and 
were sensitive from 20 keV - 1 MeV. Also included in each 
detector module was a spectral detector (SD) made of a 
smaller 127 cm^ Nal(Tl) crystal. The SDs were optimized 
for energy coverage and resolution, and were operated over 
the range 10 keV - 100 MeV. An illustration of CGRO and 
a BATSE module can be seen in Fig. 

Data from the detectors were recorded in several for- 
mats, of which two recorded by the LADs are of pri- 
mary interest for this work: the DISCLA data, consist- 
ing of 4 energy channels recorded in 1.024 second bins 
and CONT data, comprising 16 energy channels with a 
2.048 second resolution. Triangulation between LADs was 
used to locate GRBs crudely, to within a few degrees, de- 
pending on the stren gth of the flux ( Briggs et alJ 119991: 
IPendleton et alJll999t) . BATSE was a highly succesful ex- 
periment that detected almost 3000 GRBs in a mission 
which far exceeded the nominal lifetime of the experi- 
ment. The experiment was still providing good quality 
data when concerns over degradation in the performance 
of the CGRO attitude control system lead to the satellite 
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Fig. 2. A schematic illustrating the Earth Occultation 
technique for detecting astronomical gamma-ray sources 
with BATSE. During the CGRO spacecraft's orbit, mea- 
surable steps are seen in the BATSE detector count rates 
as a gamma-ray emitting object sets below and rises above 
the Earth's horizon. The occultation angle, (3, is defined 
as the angle between the orbital plane and the object mea- 
sured at the geocentre. 

being de-orbited safely over the Pacific in June 2000. The 
identification of GRBs with cosmological rather than lo- 
cal events can be lar gely attributed to th e wealth of data 
collected by BATSE l|Meeefan et alJUool . 

3. Earth Occultation Imaging with BATSE 

It was realised before the CGRO launch that the angular 
resolution for non-transient gamma-ray sources could be 
greatly improved, over that obtainable for GRBs, with 
knowledge of the position on the sky of the Earth's limb 
llFishman et al.lll983 IPaciesas et al.lll985HFishman et all 

A short description of Earth occultation analysis 
of persistant gamma-ray source s is given below, but the 
reader is directed to the work of lHarmon et al.l l|2002|) for 
a full discussion of the technique. 

The Earth Occultation technique is a way of perform- 
ing imaging, with simple non-pixellated or low spatial res- 
olution detectors, by using the Earth as a temporal mod- 
ulator of the gamma-ray flux. The tec hnique is possible 
with t he BATSE SDs (see, for example. iMcNamara et all 

however, due to their larger size, and hence sen- 
sitivity, data from the LADs only has been used in this 
work. The Earth, seen from a BATSE detector, has a di- 
ameter of ^--^140° and will sweep rising and setting arcs 
across the sky at the CGRO orbital period, ~90 minutes. 
As gamma-ray emitting objects set below the Earth's hori- 
zon downward steps are seen in the detector count rates 
at the orbital period (and vice-versa for rising objects, 
see Fig.|51l. The height above the background of the steps 
gives a direct measurement of the flux of the source being 
occulted. With knowledge of the detector orientation and 



the relative position in the orbit of the spacecraft with 
respect to the Earth, it is possible to project an image of 
the Earth's limb on the sky at the time of the source oc- 
cultation. An image of the gamma-ray emission from the 
sky can then be constructed by building up the images 
of the Earth's limb from many occultations. The intrinsic 
angular resolution depends on the rise time of an occulta- 
tion step. For a source crossing Earth in the CGRO orbital 
plane the transit time from full transmission to full atten- 
uation, Ai, is ~8 seconds, due to the increasing absorption 
of gamma-rays in the varying column density of the atmo- 
sphere as the source decreases in altitude. As the source 
occultation happens further away from the orbital plane 
At increases as 1/ cos /?, where f3 is the occultation angle 
measured at the geocentre between the orbital plane and 
the source. At, for a source seen by a satellite with orbital 
period P translates to an angle of 360° x At/ P. This angle 
represents the optimum angular resolution, for imaging a 
source, that is obtainable with this technique. Hence, for 
the BATSE/CGRO system the ideal angular resolution is 
~ 0.5°. When f3 > 70° no occultation is seen although the 
occultation is effectively unusable when > 68° as the 
sourc e skims across the poles of the Earth ijHarmon et alJ 

iooi). 

The inclination of the CGRO orbit was 28° with a pre- 
cession period of 53 days. This orbit is such that the whole 
sky undergoes occultation during the precession period, al- 
though regions of sky near the poles with |41° < 5 < 82° | 
receive less coverage than the rest of the sky. In total, in- 
cluding effects for lost data etc., between 50% and 90% of 
the entire sky can be observed in any 53 day period. The 
Earth occultation technique, when used with an instru- 
ment such as BATSE with a large field of view is a very 
attractive tool for the production of all-sky images. 

3.1. Standard Occultation Analysis Techniques and 
their Limitations 

Earth occultation techniques have been applied to the 
BATSE database since the launch of CGRO to observe the 
gamma-ray sky. Techniques were developed by workers at 
the Marshall Space FHght Centre (MSFC, note that the 
science groups of MSFC were moved to the National Space 
Science and Technology Centre (NSSTC) in 2001) to al- 
low fast monitoring of bright transient or flaring sources 
by searching thr ough the count rates for occultation step- 
like features fsee lHarmon et aP (^2002') for a full explanan- 
tion). By measuring the occultation time the location of 
the source could be deduced, from the orientation of the 
Earth's rising and setting limbs, and a position determined 
to ^ 1° accuracy. Examples of the use of this technique 
are the discovery of t he gamma-ray sour ces GRS 1716- 
249 teallet et al.ll993|) and GRO 0422-^32 l)Harmon et alJ 
1992bir and observatio ns of the X-ray nova 4U 1543-47 
(iHarmon et alJll992a|) . However, the sensitivity of the 
step searching technique is limited by the difficulty in ac- 
curately determining the background over several orbits 
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and by non-statistical noise introduced by bright pulsat- 
ing sources with periods of orde r a few minutes, such as 
Vela X-1. iHarmon et alJ l)2002() quote a 3cr sensitivity of 
~ 200 mCrab for a source with a spectrum like that of the 
Crab Nebula (4.6 x 10^9 ergs cm'^s-i for 30 - 250 keV). 

For known sources a fitting procedure can be used, on 
a 20 second fitting window, applied within a short (220 
or 240 second) section of data, centred on the expected 
occultation of the source in question. For each LAD and 
energy channel the model fit consists of a quadratic curve 
describing the background and a set of response vectors 
that model the effects of all sources contained within the 
fit window (i.e. the target source and any other interfer- 
ing background source). This allows multiple steps to be 
combined, giving an improvement in sensitivity, over the 
step search technique, proportional to the square root of 
the number of steps used. In one day it is possible to have 
between and 32 steps, and so the the improvement in 
sensitivity is about a factor of four (since the average num- 
ber of steps used is 16). In principle increasing the size of 
the fit window reduces the statistical uncertainty in mea- 
suring the step size and hence improves the sensitivity. 
However, the length of fit window is constrained heavily 
by the ability of the background model to describe vari- 
ations in the count rate. As the fit window increases a 
more complex model may be needed to account for slower 
varying terms in the background as well as accounting 
for the higher probabilty of encountering other steps from 
faint, unknown or flaring sources. Indeed, the effects of 
interfering sources could a ccount for a 30-50% error in de- 
termi ning BATSE fluxes l|Westmorel I2OO2I: iMalizia et all 
Hence a compromise must be made to reduce sys- 
tematic uncertainties in the fit at the expense of increased 
statistical error. 

The step searching and fitting techniques provide a 
flux and a position but are not capable of producing im- 
ages. In particular the only consistent way of dealing with 
the problems of interfering sources in the flt window is 
to develop a generalised, unbiased imaging capabilty. The 
timing information about occultation steps in the count 
rate light curve can be transformed into spatial informa- 
tion for the construction of images. The time derivative of 
the count rate in a BATSE detector is equal to the line 
integral of the flux along the Earth's limb. During an oc- 
cultation the flux is sampled at discrete time intervals as 
the occulting part of the Earth's limb tracks across the 
source at a constant angle, 9{(3) (see Fig.EJ. For a single 
location on the sky, the flux will be sampled at different 
angles for subsequent rises and sets of the source (except 
in the case where /3 = 0) and as the CGRO orbit precesses, 
allowing the reconstruction of the image. The occultation 
transform imaging technique uses a combination of Radon 
transformation and maximum entropy methods to make 
imag es from BATSE occultation data l|Zhang et alJIlOQSl 
Il995|) . The technique has proved successful in locating sev- 
eral new transient sources with better t han 0.5° accuracy 
I Lapshov et al 1 ll993HZhang et a,lJll994l ) and also allowed 
the observations of weak emission from low mass X-ray 



binaries, active galactic nuceli and supernova remnants 
llBarret et al.lll996tlMalizia et alJlioOoHMcCoUough et alJ 

The occultation transform imaging technique suffers 
from a number of problems that limit its usefulness for 
observing weak sources and crowded regions, and make it 
ineffective for the production of all-sky images. Firstly the 
use of iterative non-linear reconstuction techniques, such 
as the maximum entropy method are significantly compu- 
tationally intensive. This makes the consideration of long 
sections of data impractical and it is unclear how images 
constr ucted from a num ber of short sections can be com- 
bined llWestmorell27)0l . Also, the Radon transformation 
uses an implicit assumption that the Earth's limb is lin- 
ear and so only a small part of the occulting arc can be 
considered. This effectively limits the field of view of the 
final image to 20°x20°, since larger images produce dis- 
tortion of point source locations near the edges of the field. 
The Radon transformation, as applied to the occultation 
transformation, states that an image M{x, y) can be re- 
constructed from a set of g{9^ t) corresponding to the time 
series of the parallel line integrals, inclined at an angle 9 
and swept across M{x,y) in time t. In principle, a simple 
point source can be reconstructed from two orthogonal 9. 
However, in practice, the object may not be a point source 
and there is likely to be more than one source in the field 
of view. The effects of systematic and random noise mean 
that 9 must be well sampled in order to accurately recon- 
struct the image. The sampling of t is also incomplete if 
the field of view of the reconstructed image is restricted. 
A bright source just outside of the reconstructed field of 
view will make a contribution to the background flux in 
the image. Fig. shows the projection on the sky of im- 
ages of the Earth's limb seen by BATSE during occulta- 
tions of the Crab Nebula, built up over a 53 day period. 
It can be clearly seen that a source has an influence over 
a very large part of the sky image and also shows that, 
even over a whole precession cycle of the CGRO orbit, 9 
is not completely sampled. 

4. All-Sky Imaging with Maximum Likelihood 

To fully exploit the A-k FOV of the BATSE instrument 
it is necessary to take a different approach in reconstruct- 
ing the image to the deconvolution used in the occultation 
transform method, outlined in Sec. 13. II which is limited to 
small areas of the sky. The Likelihood Imaging Method for 
BATSE Occultation Data (LIMBO) code is a collection of 
computer programmes developed at Southampton to cre- 
ate all-sky maps from the BATSE data and is based on 
initial work by J. Knodlseder of Centre d'Etudc Spatiale 
des R ayonnements (CESR) ()Knodlseder,.1999 : Shaw et alJ 
l200lh . A pixellated image is btiilt up by using a Maximum 
Likelihood Ratio (MLR) test for source emission from a 
given source location on a predefined grid of sky positions. 
The problem of image contamination by sources outside 
the image FOV, seen in the occultation transform tech- 
nique is avoided by this method as the total dataset, com- 
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Fig. 3. The projections on the sky, in galactic coordinates, of occultation limbs for sources seen by BATSE during a 24 
hour period: a) The Crab Nebula for TJD 09586; b) The Crab Nebula for TJD 09560; c) Cygnus X-1 for TJD 09586. 
The time difference between TJD 09560 and TJD 09586 is approximately 50% of the CGRO orbital precession period. 
Note how the quality of positional information in two dimensions varies with source position and time. The result of 
adding limbs for a full 53 day precession cycle is shown in d) and shows how the Crab Nebula position is well defined 
with no positional ambiguity caused by secondary intersections of the limbs, such as is seen in a) and c). However, it 
is clear that the angular sampling of the source is still not complete. The crosses mark the positions of the sources. 



prising the whole sky image, is used in the calculation of 
the MLR for each pixel. It should be stressed that this 
technique differs from the occultation transform method 
in that it produces maps of the results of a number of sta- 
tistical tests for emission rather than a deconvolution of 
the data. There are two main parts to the LIMBO imag- 
ing process, the background corrected data is first trans- 
formed through a differential filter before the likelihood 
tests are applied and the maps created. 

4.1. Flat Fielded Data 

When considering long sections of BATSE data a good 
understanding of the highly variable gamma-ray back- 
ground encountered by the instrument in orbit is essen- 
tial to enable the reduction of systematic er rors intro- 
duced in flux determination. iLing et aD l)2000f) of the Jet 
Propulsion Laboratory (JPL) group report the use of es- 
sentially similar analysis techniques to MSEC, but with 
a multi-parameter background fit based on a combina- 
tion of empirical and semi-empirical functions derived, 
in part, from data taken by the spacecraft. This allows 



them to extend the leng t h of d ata to be fitted to 1 day. 
However, iHarmon et al.l l|2002|) note that there may be 
limiting systematic uncertainties remaining, which man- 
ifest themselves as significant uncharacteristically hard 
emission being reported in a number of the published JPL 
source spectra. 

The BATSE Mass Model (BAMM) is a Monte- 
Carlo based background model and was developed at 
Southampton primarily for use with the LIMBO tech- 
nique. Information about the BATSE detector orienta- 
tions, the position of CGRO in orbit in relation to the 
Earth and the local magnetic and radiation environ- 
ment are taken from the data stream and used as the 
basis for a physical background calculation based on a 
database of simulated particle fluxes. By considering the 
energy dependent contributions to the background from 
three components (cosmic diffuse and atmospheric albedo 
flux gamma-rays, prompt cosmic ray interactions). The 
BAMM model is capable of reducing the variations in the 
BATSE background flux, by a factor of 8-10, simply by 
subtracting the modelled flux from the data. This effec- 
tively nullifies all systematic effects on the imaging anal- 
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Fig. 4. The BATSE CONT LAD 1 25 - 35 keV count rate for TJD 9586 (shown in grey) with corresponding results 
from the BAMM model overlaid (dark curve). BAMM reproduces the shape, amplitude and variations in the BATSE 
count rates very accurately: No further fitting or renormalising of the model has been performed in this graph, except 
for the addition of a constant offset of 351.5 counts per bin to the model for display purposes only. The counts are 
shown per 2.048 second CONT time bin. 
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Fig. 5. Top: The result of passing the flat-fielded BATSE CONT LAD 1 25 - 35 keV count rate for TJD 09586 through 
the differential filter of Eq.^ The counts show clear peaks, due to the occultations of bright sources. The peak positions 
can be compared with the bottom graph, which shows the expected differential response vector multiplied by the value 
of the flux scaling parameter, a, for a source at the Crab Nebula position. The gaps and shortening of some response 
peaks are due to missing data in the CONT stream. In both graphs the filter applied was fa — 4:, fi, — 14. 



ysis, producing a 'flat-fielded' data set. A full description 
of BAMM is outside the s cope of this document but is 
the s ubject of other works ijShaw et alJl2003t iDean et aD 
l200nl) . A general discussion of mass modelling techniques 
for appli cations in h i gh en ergy astrophysics can also be 
found in lDean et alJ ((2003f l. An example of the effective- 
ness of BAMM in accurately modelling the BATSE count 
rates can be seen in Fig. 01 



4.2. Differential Filter 



As CGRO orbits the Earth a path of sky of width up to 
140° may be occulted by either the rising or setting limb. 
Hence it is likely that many sources of varying strengths 
will be observed, leading to a superposition of steps in 
the count rate. The useful information about the source 
flux is contained in the short section of data surrounding 
the step itself, therefore it is attractive to deal not with 
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the absolute detector counts but with their time rate of 
change. 

The flat-fielded BATSE data is passed to a simple dif- 
ferential filter, given by 



E3=i+lb 
i=i 



i+fa 



h-fa + 1 



(1) 



The filtered data vector, o, is transformed from a count 
rate, r, containing a superposition of occultation steps to 
a set of occultation peaks by the filter with inner and 
outer bounds given by fa and /& respectively (see Fig.|Sl). 
A square step at bin i will be transfromed to a triangular 
peak of width 2/, where f = ft — fa- The choice of / has 
implications for the sensitivity and angular resolution of 
the imaging technique; larger / will give a more signifi- 
cant occultation peak and improve sensitivity, smaller / 
returns less significant but sharper peaks that will lead 
to better angular resolution. The filter properties are dis- 
cussed further in Sect. 14.61 

The raw BATSE data contains occasional gaps of null 
data where, for example, the instrument is turned off dur- 
ing passage through the South Atlantic Anomaly (SAA) 
or where contact is lost between CGRO and its telemetry 
relay satellites. Left untreated the gaps cause large artifi- 
cial peaks in the filtered data near each gap, which mimic 
sources. The result of differentiating bin i just before a 
gap is a large positive peak since Eq.^ depends on counts 
from bins i + f, at least some of which will be null data 
from within the gap. One way to reduce this effect would 
be to ignore all data within / bins of a gap, although this 
may remove ~ 10% from the total of good data available. 
Instead a simple linear interpolation is used to pad the 
gaps in the flat-fielded data, this stops the filter blowing 
up in gap-data boundaries and conserves the maximum 
of available data. The padded data from within the gap is 
marked as bad in the filtered data set and not used further 
in the analysis. 

4.3. Maximum Likelihood Imaging 

A MLR test is used to determine the significance of flux 
detection at a single point in a given grid of source po- 
sitions from the filtered data vector a. The MLR, A, is 
calculated from the difference 



A — Co Cs 



(2) 



where Co and Csrc are the statistics for the null 
(no source is fitted to o) and source (a source with free 
strength is fitted to o) hypotheses respectively: 



Co = 



aei) 



(3) 
(4) 



The weighting factor, di, is taken as the statistical er- 
ror on Oi. This becomes the average poisson error on the 
filtered BATSE count rate. Assuming that the error on 
the modelled background is small; 



2/ + 1 



(5) 



The response vector, e, describes the expected posi- 
tion, shape and amplitude of the peaks as a function of 
source strength and position on the sky. It is also differ- 
entiated so that the series of expected steps in the count 
rate for each source position are transformed into peaks, as 
shown in Fig.|Sl The differentiated response scales between 
(no change in transmission) to ±1 (maximum change in 
transmission) as a function of the angular response of each 
LAD. The scaling factor, a, is proportional to the flux re- 
ceived from each point in the sky for all eight LADs, 



(6) 




and has a statistical uncertainty given by, 

/ \ -1/2 



p) 



The distribution of A is the same as xt with v degrees 
of freedom. When searching for a known source v = I and 
hence the detection signiflcance, in gaussian a, is given 
by -^A. Because of this relationship gives source sig- 
nificance directly. Since a satisfies the condition of min- 
imum it can be considered as a property of the data 
and the response vector and it is not a free fit parameter. 
The test is easily extended to searches for new sources of 
gamma-rays at unknown positions; two extra degrees of 
freedom are added for the {x, y) sky position and so = 3 
l Knodlsedeillll99<j> ). 

The likelihood image, A(x,?/), is built up from the su- 
perposition of the arcs caused by the occultation of point 
sources by the Earth's limbs. The effect of the occulta- 
tion angle, (3, on the angular resolution of the imaging 
technique was discussed in Sect. |21 A second effect im- 
proves the resolution at moderate [3 when more than one 
occultation arcs are considered. At (3 ~45° the rising and 
setting limbs are almost orthogonal, and give good posi- 
tional information in both the x and y directions of the 
map simultaneously. In addition the natural precession of 
the CGRO orbit means that a number of angles are sam- 
pled every 53 days. 



4.4. LAD Response Matrix 

The calculation of response matrices to describe the sen- 
sitivity of gamma-ray detectors as photon energy and in- 
cident angle changes is a difficult task. The app roach used 
here is based on work bv lHarmon et alJ ( 20021) . The stan- 
dard occulation technique was used to observe the Crab 
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Fig. 6. The flux sensitivity (3(t) and angular resolution 
(FWHM) of the LIMBO technique for varying filter sizes, 
/, measured in number of 2.048 second CONT bins. The 
measurements are based on one days observation of the 
Crab Nebula imaged with the LIMBO technique over a 
0.1° grid of positions. 

nebula over a considerable portion of the CGRO lifetime. 
During this time the Crab was seen by all eight detectors 
at various angles incident to their optic axes, 9. A func- 
tion was then fitted to the points for each LAD and energy 
channel assuming azimuthal symmetry of the instrumen- 
tal response. The response, R, is described by an equation 
containing 4 parameters, Ai, A2, A3 and A4, which are 
unique for each energy channel and detector. The eq ua- 
tion for the response, taken from lHarmon et aD l)2002j) . is 
as follows: 



Al cos(6')exp'' 



os(8) sin(6l + A4 



(8) 



This calibration gives maps of a in units of the nor- 
malised Crab Nebula response. 

4.5. Combination of Images 

The sensitivity of the imaging technique can be increased, 
by approximately the square root of exposure time, by 
considering longer observations. It is possible to run the 
maximum likelihood imaging technique, of Sect. 14.31 on 
any length of time required. However, in practice, the 
maximum length of time that can be considered in this 



way is limited to a few tens of days by the availability of 
computing resources, since large amounts of data need to 
be read into memory at once. Also, once the image has 
been made, it is not possible to simply extract an image 
from a subset of the data without running the analysis on 
that data alone. Hence it is attractive to be able to make 
images of short periods of time (e.g. daily) and combine 
them into longer, more sensitive observations at a later 
date. This would also allow the time variability of bright 
sources, found in the high sensitivity image, to be investi- 
gated on shorter timescales. 

It is not possible to add individually produced likeli- 
hood images together directly. This is seen from the de- 
pendence of Eq.01on Eq.El since a is the ratio of two sums, 
formed from the complete data set in question, namely; 



D 



(9) 



(10) 



It is possible, however, to calculate the total flux value, 
ax, for each image pixel by summing the values of F and 
D from shorter time periods. In this case d will be used 
to index a number of images, each made from single, day 
long, observations. It is important to emphasise, since even 
for a constant source of gamma-rays, the values of may 
be different from image to image (i.e. due to data gaps, 
changing orientation of the spacecraft, etc.), that 



ax 




ad 



(11) 



Hence, if a data set, indexed by k = {1, . . . ,d,d + 
1, . . . , N}, is split into two sets, with Csrca ^^'^ Q^a eval- 
uated for k = {1, . . . , d}, and Cgrcb and at evaluated for 
k — {d + I, . . . , N}, then it is clear that Csrca + Csrc^ 7^ 
CsrcT because 



i(k=d) 

E 

i=i(k=i) 



{Oi - aaCiY 



{oi - abCiY 



+ E 

i{N) 2 

[Oj - arej ) 

i=i{k=\) * 



(12) 



Images of statistical significance of gamma-ray emis- 
sion can be made by scaling (Xt by its statistical uncer- 
tainty, Aa, given by Eq. [3 With this method it is also 
possible to combine images made from individual BATSE 
energy channels, which allows the sensitivity of the imag- 
ing technique to be increased whilst retaining spectral, as 
well as temporal, resolution. 

4.6. Differential Filter Choice 

The choice of differential filter width, /, has implications 
for the sensitivity and angular resolution of the imaging 
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technique. Choosing a larger / will give a more significant 
occultation peak and improve sensitivity, smaller / returns 
less significant but sharper peaks that will lead to better 
angular resolution. Fig. shows the performance of the 
imaging technique for various choices of /. The figure is 
based on a 53 day observation of the Crab Nebula using 
one BATSE energy channel (~ 25 — 35 keV) and an image 
grid resolution of 0.1° . The sensitivity improves roughly 
as l/V/ but is seen to level out as / increases to the 
point where the timescale of variations in the background 
approach / x 2.048 seconds. For all choices of / the peak 
flux in the map was located within < 12' of the actual 
Crab position of = (184.56° ,-5.78° ). It can be 

seen that / w 10 — 20 offers a good compromise between 
sensitivity and resolution. It is clear however that different 
/ may be used to optimise the technique for particular 
science goals. 

It is possible to make an improvement to the sensi- 
tivity of a particular chosen / whilst retaining reasonable 
angular resolution. The size of a differentiated occulation 
peak depends on the difference between the mean counts 
calculated before and after the time bin in question. Since 
the occultation steps are not infinitesimal (they take about 
10 seconds), the mean counts contain a contribution from 
the transition of the source from full transmission to full 
occultation that acts to reduce the average difference be- 
tween the two means. This causes a reduction in sensitiv- 
ity which becomes more important as / decreases towards 
the step time of about 5-6 bins. This can be avoided by 
adding an inner bound to the differential filter. By making 
a small gap in the centre of the filter systematic errors in 
calculating the mean counts either side of a step are re- 
duced but angular resolution can be retained. Fig. [3 and 
TableQlshows the effect of applying several different filters 
to a simple model of an occultation step. Clearly / = 10 
with a small inner gap offers the best combination of an- 
gular resolution (given by the HWHM of the curves) and 
sensitivity. 

Table 1. Approximate sensitivity and angular resolution 
of differential filters used in Fig. |7| 



/ 


/. 


Sensitivity 


Angular resolution 


(bins) 


(bins) 


(%) 


(°) 


10 





80 


0.75 


10 


2 


97 


0.89 


10 


5 


100 


1.21 


20 





90 


1.38 



A possible extension to the differential filtering tech- 
nique may be to use different resolution maps to build up 
images of extended emission. By choosing different val- 
ues of / that optimise the coverage of spatial frequencies 
in the image it may be possible to investigate, for exam- 
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Fig. 7. A simple model to show the effect of various differ- 
ential filters on an occultation step of duration 12 seconds 
(solid line). Filters with / = 10 were applied with inner 
gaps, fa, of (dotted), 2 (short dash) and 5 (long dash) as 
well as a filter with / = 20 and inner gap = (dot-dash). 

pie, diffuse gamm a-ray emission from the Galactic Plane 
l|Westmorell200^ . 

4.7. Maps of statistical significance of source emission 

An all-sky map for 489 days in the period TJD 09448 
- 09936, produced from 7 CONT energy channels (25 - 
160 keV), is shown in Fig. [S] Using simple statistical ar- 
guments, and neglecting any systematic background, the 
map has a 3cr sensitivity of 3 mCrab. This can be consid- 
ered as the first all-sky gamma-ray image to be produced 
since the HEA01-A4 mission. Although it is easy to see 
the brightest sources in the sky, such as the Crab Nebula 
and Cygnus X-1 it is clear that the imaging technique 
is heavily affected by artefacts stemming from the bright 
sources. This is a particular problem in the crowded re- 
gion of the Galactic Centre, where source confusion is a 
problem. 

The artefacts seen in the figure are images of the 
Earth's limb from individual occultations. An individual 
limb implies a probability of a source being observed any- 
where along its length. The source location accuracy im- 
proves as more limbs are superimposed and the size of the 
likelihood peak increases. Although the brighter sources 
can be located accurately from their likelihood peaks the 
artefacts from limbs passing through the source position 
may be more significant than the likelihood peaks found 
at the position of fainter sources. 

5. Image Artefact Removal I: Removal of 
Individual Sources 

As mentioned earlier the diameter of the Earth seen from 
CGRO is about 140°. Hence the limbs crossing through a 
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Fig. 8. Map of the significance of 25 - 160 keV flux from the whole sky for a total of 489 days from TJD 09448 
to TJD 09936 made using LIMBO. The colour scale of the map has been limited to a quarter of the peak in order 
highlight details at lower significance. Many well known, bright sources are visible and the peak significance of 920 a 
is consistent with the position of the Crab Nebula. There is a limited systematic across the whole image but there are 
serious artefacts which limit the usefulness of the image for observing weak sources or crowded regions. 



source position can have an effect over a large region of 
the sky. This is especially obvious from the map of Fig. |S1 
for the bright sources, such as the Crab and Cygnus X- 
1. The effects of these sources can be removed from the 
data by making two passes of the likelihood calculation. 
By first calculating the response function for a particular 
source position, the effect on the entire image of a source 
at this position can then be removed by adding the inverse 
of the response to the likelihood calculation. Fig. shows 
the removal of Cygnus X-1 and all its associated artefacts 
from an image. When the confusing influence of Cygnus 
X-1 has been removed it is possible to see a faint source 
that is consistent with the position of Cygnus X-3. 

The iterative removal of individual sources in this man- 
ner may not, however, be an ideal method of reducing the 
influence of image artefacts so that fainter sources may be 
observed: The response measured from each source posi- 
tion will include a systematic effect due to artefacts from 
sources at other positions. This will lead to a cumula- 
tive error in estimating the flux each time a source is 



removed, although the seriousness of this effect may be 
limited if only the most significant of sources are removed 
first. Finally repeated re-running of the LIMBO code will 
clearly incurr overheads in computing time that may be 
prohibitive. 

6. Image Artefact Removal II: Image Cleaning 
Algorithm 

All pictures or maps produced by any imaging technique 
consist of the true spatial distribution of light (the object 
function) convolved with some distorting function specific 
to the instrument or technique being used (the impulse 
response). In optical astronomy, the impulse response is 
minimised by carefully designing the instrument to re- 
duce aberrations leading to a compact, symmetrical Point 
Spread Function (PSF). At other wavelengths, where in- 
strument designers have less control over the optical de- 
sign, a post observation approach must be applied to re- 
duce distortions in the image. One such example is in ra- 
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Fig. 9. An example of the removal of source artefacts by making a second pass of the data with the likelihood function 
adjusted to include an inverse of the response at the position of Cygnus X-1 (maked by 'x'). The left image shows the 
LIMBO all-sky map for data from TJD 9450 - 9900, 25 - 35 keV. In the right image Cygnus X-1 has been removed. 
Note the faint response near Cygnus X-1 that is apparent in the right image, which is consistent with the position of 
Cygnus X-3 (marked by '-I-'). 



dio in terferometry where the CLEAN algorithm l|Hogboml 
I1974I) has been used for many years to remove artefacts 
from radio maps introduced due to the finite number of 
radio dishes, the distribution of the baselines used and the 
motion of the Earth relative to the source l , Tan 1986i() . 

The CLEAN algorithm works by defining two data 
spaces; the raw map Mn (consisting initially of the original 
sky image) and the clean map Mc (of the same dimen- 
sions as the raw map, but initially empty). The algorithm 
proceeds by first finding the brightest pixel in Mr and 
generating the expected PSF for this position. A small 
fraction of the PSF, defined by a gain set by the user, is 
removed from Mr. The total flux that has been removed 
is integrated over the PSF. If required, this flux can than 
be convolved with an ideal PSF, such as a Gaussian dis- 
tribution with a realistic width, and added to Mc at a 
position centred on the position of the brightest pixel in 
Mr. This process is repeated iteratively until some pre- 
defined criterea, such as the minimum significance of the 
brightest pixel remaining in Mr, is reached. At this point 
what remains of Mr should resemble intrinsic noise, which 
can be added to Mc to produce a more realistic final im- 
age. 

At any time the shape, size and orientation of a 
BATSE occultation limb can be calculated given informa- 
tion about the orbit of CGRO and the detector pointings. 
Since the artefacts produced by point sources in LIMBO 
are due to the accumulation of occultation limbs, an im- 
pulse response function may be calculated for any posi- 
tion in the sky based on the knowledge of the history of 
the position of the Earth's limb. This response to a point 
source is conceptually similar to a PSF for the observa- 
tions, although it should be highlighted that the physical 



interpretation differs from the PSF used in other imaging 
modalities. For example, in optical astronomy the PSF 
represents the characteristic way the imaging system re- 
distributes the point source of flux in the image; a point 
source of photons in the sky becomes extended in the im- 
age. In LIMBO the PSF represents the characteristic way 
the imaging system redistributes the probability of the es- 
timated position of the source; a point source in the sky 
produces an image that is consistent with emission from 
a number of different, but mutually exclusive, positions. 

The PSF for any particular sky position is highly com- 
plex, as is shown in the example in Fig.O However, if the 
PSF of a selected source position can be calculated it is 
possible to subtract that source's contribution to the total 
image. 

6.1. Calculation of the UMBO Point Spread Function 

The PSF for the LIMBO technique is highly complex as 
it is defined by the set of loci of the Earth's limb for any 
given source position and depends on the orientation of 
the detectors with respect to the position of the spacecraft 
in the orbital precession cycle. Consequently the PSF for 
any given source is temporally, spatially and spectrally 
variant. It is possible to make an accurate determination 
of the PSF by creating a dataset from the known response 
vector of a target source and applying the LIMBO code to 
it. However, this is an extremely inefficient approach as in 
an all-sky image the majority of pixels for which LIMBO 
generates a response and evalutes the MLR never undergo 
an occultation simultaneously with the target source. A 
more efficient approach has been developed bv 'West morel 
((2002) , which uses knowledge of the position of the Earth's 
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limb over time to restrict the number of pixels LIMBO 
tests for a response from. For each time bin during the 
occultation, the pixels through which the Earth's limb 
passes can be identified by using the spherical polar coor- 
dinates of the Earth's centre (as seen from CGRO), which 
are recorded along with the counts, in the BATSE data 
stream. Once all of the relevant pixels have been identified 
LIMBO can be run on this smaller set. For a 2° resolution 
grid of sky positions this approach leads to a reduction in 
the number of pixels for which the MLR is calculated from 
16,471 to ~ 1000, which allows a significant reduction in 
the amount of computing time required. 

6.2. Choosing a Linear Cleaning Parameter 

An essential requirement for the CLEAN algorithm to 
work is that the imaging property being reduced is linear, 
i.e. the quality of the final map should not be a function of 
the order in which contributions to the background from 
individual sources were subtracted. This is not true for 
the MLR, which is based on a geometrical sum. However, 
F, (recall Eq. can be shown to be linear for an image 
containing a superposition of S sources. The differential 
filter of Eq. ^can be written in a simplified form as 



1 



(13) 



If the filter is associated with an operator, O, such that 



0(rj) = Oi, and rj = ^sj is substituted into Eg. 1131 

it can be shown, since r is linear, that 



/ + 



^EE^^. = Eo(' 



(14) 



Hence the law of superposition, 0(J2 



s=S 



^^"-^ 0{rsj), is true and it has been shown that the dif- 
ferential filter is a linear operator. Since both and 
are constants it follows that F is also linear, 



F 



EE 



0{rsj)ei 



(15) 



It does not follow that the flux maps are linear. 
However, since D (Eq. I10|l . is a constant it is possible 
to run the cleaning algorithm on F and use the results of 
this to derive the cleaned flux and log- likelihood maps. 

6.3. Identifying which Pixel to Clean 

For some systems the identification of the sky pixel loca- 
tion at which to clean, P(i^b)i is trivial - it is simply the 
brightest pixel. For the LIMBO system this is not neces- 
sarily the case since the peak in F is not necessarily in a 
pixel consistent with the true source position leading to 
calculating the PSF for an incorrect position. For this rea- 
son P(i^b) is found by first generating the PSF for the most 
significant pixel and then cross-correlating this on a grid 



surrounding the most significant pixel. The pixel position 
with the highest correlation coefficient is then assumed 
to be the correct pixel. A second PSF is regenerated for 
this position and a small fraction of it subtracted from the 
image and added to the clean map. 

Choosing the incorrect P[ifi) is particularly problem- 
atic for the LIMBO images since, for these initial tests of 
the technique, the pixels are 2°x 2°and hence the error 
can be great. Even if the correct pixel is chosen the PSF 
will still not be generated at exactly the correct position 
since it must be generated at the centre of the pixel; the 
maximum discrepancy being + 2^ = 2.8°. In order to 
address this problem two strategies are introduced. Firstly 
there is a pre-cleaning stage where the majority of the 
contributions from the brightest known sources, with well 
defined positions, are removed. Secondly, a catalogue of 
probable sources is checked as each p^i^b) is found. The pre- 
cleaning stage is particularly important since the system- 
atic error introduced into the final Mq map by cleaning 
with the incorrect PSF is dependent on the total amount 
of flux removed with that PSF. Hence, if bright sources, 
such as the Crab, are cleaned with only a slightly incor- 
rect position, the resulting errors will be large compared 
to the faint sources whose detection is the ultimate goal. 

For a list of known sources, the PSF at the exact po- 
sition each of the sources is generated and used prior to 
running the cleaning algorithm to remove most (typically 
95%) of the contribution to Mr. As new sources are dis- 
covered they may be added to this list and the cleaning 
algorithm run again to find the next level of sources. As 
the cleaning algorithm progresses each P(i^b) CLEAN pixel 
is also compared to a list of known sources that it is be- 
lieved, though not definitively known, to be hard X-ray 
sources. This list comprises, amongst others, the HE AO 
catalogues, known AGN, and transient Galactic sources. 
If within P(i^b) there is a candidate source then the PSF 
may be generated at the exact position of this source. Care 
must be taken however in crowded regions. 

7. LIMBO Image Cleaning - Preliminary Results 

The results of applying the cleaning algorithm described 
above can be seen by comparison of Fig. ^| and Fig. |S1 
Sources were extracted from the map us ing the SExtractor 
2.2.2 (SE) software, developed by 'Bcrtin fc Arn outi 
with the following settings: DETECT.THRESH 
set to 1 sigma; CLEAN_PARAM (cleaning efficency) of 
1.0; a wavelet pass-band filter, mexhat_l.5_5x5.conv. The 
fluxes and their associated errors are measured using SE's 
ISO_CORR aperture option. 

Two maps were produced to use with SE; one con- 
taining the cleaned components only and one with the 
residual flux, remaining after cleaning, added back to the 
cleaned components. The latter map is of course the more 
scientifically realistic map and when running the SE code 
directly on it produced a catalogue of detections in the 
image, shown in Tabled The fiuxes of the sources in the 
table have been normalised to the Crab flux. 
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Fig. 10. The results of performing the image cleaning technique on the all-sky map of Fig. |S| See Tables |21 and |31 for 

a summary of the identified sources from the image. 



The sources extracted by SE were compared against 
the INTEGRAL Scie nce Data Centr e (ISDC) catalogue of 
high energy sources ijEbisawa et al.li2003afij) . which con- 
tains over 1000 sources with flux estimates in the 20 - 
200 keV band. Points of emission in the map have been 
identified with sources from the catalogue, based on their 
positions and estimated brightnesses. It should be noted 
that, due to th e inhere nt variability of gamma-ray sources, 
lEbisawa et alJ 1 2003a|^ report large uncertainties in some 
of their flux estimates. In addition, the catalogue and 
LIMBO spectral bands are not equal and so, in some cases, 
these source identifications are subjective and should be 
treated with caution. Any source which lies within 2.8° of 
its catalogue position is within the diagonal size of one im- 
age pixel. The angular resolution of the LIMBO images is 
hmited by the 2° grid spacing. Hence there are some cases 
where it is not possible to associate the detection with 
only one known source. Source confusion is always likely 
to be a problem, especially in the Galactic Centre, with 
a technique such as Earth Occultation Imaging. However, 



recall from Sect. O that it is possible to improve the angu- 
lar resolution of the LIMBO technique to 0.5°. 

A number of sources visible in Fig. 1101 were not, how- 
ever, detected by SE. Table O shows a list of sources not 
included in Table |21 which was compiled by using the po- 
sitions where PSFs were generated during the cleaning 
process as a guide for SE to perform its source extrac- 
tion from the cleaned -I- residuals flux map. This method 
assumes that the cleaned detections are real. It should 
be noted that fake sources, added at random positions 
in the input catalogue, will return essentially zero flux. 
Hence guiding the source extraction in this manner does 
not produce false detections. The measured fluxes of these 
sources indicate the possibility of sources at that position, 
but as they cannot be independently measured may not 
be entirely trusted. It should be noted that the positional 
error for many of these source identifications is larger than 
for those in Tabled In the Cygnus region this analysis re- 
sults in the detection of EXO 2030-375 in addition to Cyg 
X-1 and Cyg X-3. 
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Table 2. List of sources detected by BATSE with LIMBO in the 25 - 160 keV band for the period spanning TJD 
09448 - 09936. The fluxes have been normahsed to the Crab Nebula flux. 



Measured Position Source Name Known Position Position Error BATSE Flux 



(l,b) 




(l,b) 


(degs) 


(mCrab) 


(71.9, 3.8) 


Cygnus X-1 


(71.34, 3.07) 


0.9 


1142 ± 18 


(-175.9, -5.9) 


Crab 


(-175.44, -5.78) 


0.5 


1000 ± 17 


(0.0, 6.9) 


GRO J1719-24 


(0.18, 7.02) 


0.2 


214 ± 9 




H 1705-250 


(-1.41, 9.06) 


2.6 




(-12.0, 2.0) 


4U 1700-377 


(-12.2, 2.2) 


0.3 


209 ± 8 




GX 349+2 


(-10.90, 2.74) 


1.3 




(-2.0, 0.9) 


GRS 1734-292 


(-1.11, -1.41) 


1.0 


201 ± 9 




H 1743-322 


(-2.87, -1.61) 


2.7 




(0.0, 0.2) 


XTE J1748-288 


(0.68, -0.22) 


0.8 


180 ± 8 




IE 1740.7-2942 


(-0.88, -0.11) 


0.9 






GX 359+2 


(-0.43, 1.56) 


1.4 






GX 3+1 


(2.29, 0.79) 


2.4 




(-16.0, 2.0) 


OAO 1657-415 


(-15.64, 0.31) 


1.7 


128 ± 7 


(0.1, 24.0) 


Sco X-1 


(-0.90, 23.78) 


0.9 


114 ± 7 


(-21.5, -4.7) 


GX 339-4 


(-21.06, -4.33) 


0.6 


93 + 6 


(-96.3, 3.7) 


Vela X-1 


(-96.94, 3.93) 


0.7 


83 + 6 


(16.7, 1.4) 


GX 17+2 


(16.44, 1.28) 


0.3 


77 + 8 


(122.0, -30.0) 


Mkn 348 


( 122.27, -30.91) 


0.9 


51 + 8 


(4.6, -1.4) 


GX 5-1 


(5.07, -1.02) 


0.6 


46 + 5 


(44.3, 0.0) 


GRS 1915+105 


(45.40, -0.23) 


1.1 


43 + 5 


(-2.00, 6.0) 


GRO J1719-24 


(0.14, 6.99) 


2.4 


42 + 5 


(-60.00, 0.0) 


GX 301-2 


(-59.90, -0.04) 


0.1 


42 + 4 


(80.00, 0.6) 


Cygnus X-3 


(79.84, 0.69) 


0.2 


36 + 5 


(-50.00, 19.5) 


Centaurus A 


(-50.48, 19.42) 


0.5 


33 + 5 


(-158.83, -4.0) 


H 0614+091 


(-159.11, -3.38) 


0.7 


31 + 8 


(-24.00, 0.0) 


H 1624-490 


(-25.08, -0.26) 


1.1 


17 + 3 


(156.00, 76.0) 


NGC 4151 


(155.08, 75.06) 


1.0 


16 + 3 


(138.00, 42.0) 


M 82 


(141.41, 40.57) 


2.9 


14 + 3 


(-70.00, 66.0) 


3C 273 


(-70.05, 64.36) 


1.6 


12 + 3 



Table 3. List of additional sources visible in Fig. ^| but not included in Table |2 These sources were extracted using 
the locations where PSFs were generated in the LIMBO analysis as an input to the SExtractor package. 



Measured Position 
(l,b) 


Source Name 


Known Position 
(l,b) 


Position Error 
(degs) 


BATSE Flux 
(mCrab) 


(-79.0, -59.0) 


ESQ 198-24 


(-88.36, -57.95) 


5.0 


62 + 8 


(-5.0, -70.0) 


NGC 7582 


(-11.92, -65.70) 


5.1 


55 + 8 


(146.1, 28.0) 


Mkn 78 


(151.10, 29.78) 


5.0 


40 + 6 


(-68.0, -38.0) 


4U 0357-74 


(-71.42, -37.29) 


2.8 


38 + 6 


(-47.9, -2.0) 


Circinus Galaxy 


(-48.67, -3.81) 


2.0 


28 + 5 


(76.8, -2.0) 


EXO 2030+375 


(77.15, -1.24) 


0.8 


22 + 5 



8. Discussion and Future Work 

The first true all-sky image of the gamma-ray sky since 
HEA01-A4 has been produced, in a single 25 - 160 keV 
energy band, using the LIMBO maximum likelihood imag- 



ing technique with BATSE occultation data. In order to 
produce the all-sky map it has been necessary to de- 
velop several new, novel techniques, such as: An accurate 
background model for the BATSE experiment (BAMM; 
■Shaw et al. (2003) ): Methods for producing all-sky images 
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from BATSE occultion data; A method for cleaning arte- 
facts from maximum likelihood images. The all-sky map 
is based on ^ 500 days of data and shows the location 
of about 25 known gamma-ray sources, including several 
extragalactic objects at 4 — Sct level (~ 15 mCrab). 

In almost all cases, LIMBO has identified emission 
from within one pixel of a source from the ISDC high en- 
ergy catalogue. However, whether all the detections from 
the LIMBO map are real soiirces of gamma-ray emission, 
that have been associated with the correct object from the 
catalogue, can not be determined at present. 

It is clear that the images presented here suffer from 
the coarse 2° resolution of the imaging grid, which makes 
identification of sources in crowded regions, such as 
the Galactic Centre, very difficult. The performance of 
LIMBO in these cases can be improved significantly by 
simply increasing the resolution of the imaging grid to- 
wards the theoretical imaging limit of 0.5°, although this 
would be at the expense of an increase in computing time. 

Since a single energy band has been considered, no 
measurements have been made of the spectra of any of the 
detected sources and fiuxes are given only in comparison to 
the detected Crab Nebula flux. The all-sky map was made 
by summing the images from 7 CONT energy channels. 
Tests have shown that there is significant imaging response 
in the individual energy channels and further information 
could be gained from summing energy channels above 160 
keV. 

It is envisaged that further investigation will be made 
of the spectral performance of LIMBO with the aim of 
calibrating the response so that the continuum spectra of 
gamma-ray sources may be examined. A rough estimate 
of the statistical sensitivity of the LIMBO technique to 
persisitent gamma-ray sources can be found by scaling 
the detection significance of the Crab nebula (990cr) in 
the 500 day image. This gives a ha sensitivity in the 25- 
160 keV range of ~5 mCrab. If the whole BATSE CONT 
dataset of ^3300 days is considered then the sensitivity is 
improved by the square root of the ratio of the observation 
times, leading to a sensitivity of ^2 mCrab. Depending on 
the brightness of the source, it should also be possible to 
make time resolved spectra for the duration of the BATSE 
mission. Assuming an equal contribution to the integrated 
flux from each of the 7 CONT energy channels and using 
the same simple statistical arguments as above, it should 
be possible to make broad spectra for gamma-ray sources 
with a constant flux of ^6 mCrab. Work continues to pro- 
duce a publically accessible database of daily all-sky maps 
at 0.5° resolution for several energy channels of the full 9 
years of available BATSE data. 
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